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Abstract 



Wc investigate level-set type methods for solving ill-posed problems, with the assumption that the 
solutions are piecewise, but not necessarily constant functions with unknown level sets and unknown 
level values. 

In order to get a stable approximate solutions of the inverse problem a Tikhonov-type rcgularization 
, approach coupled with a level set framework is proposed. Wc prove existence of generalized minimizers 

for the Tikhonov functional. Moreover, wc prove convergence and stability for regularized solutions 
with respect to the noise level, characterizing the level-set approach as a regularization method for 
inverse problems. 

We also showing the applicability of the proposed level set method in some interesting inverse 
problems arising in elliptic PDE models. 

. Keywords: Level Set Methods, Regularization, Ill-Posed Problems, Piecewise Non-Constant Coeffi- 
1 cients. 

(n: 

! 1 Introduction 

Since the seminal paper of Santosa [30], level set techniques have been successful developed and has 
\ recently become a standard technique for solving inverse problems with interfaces. As references, see 
r-H ■ [35, 17, 36, 29, 27, 8, 11, 18, 5] just for cite a few. 

. In many applications, interfaces represent interesting physical parameters (inhomogeneities, heat 

conductivity between materials with different heat capacity, interface diffusion problems) across which one 
1 or more of these physical parameters change value in a discontinuous manner. The interfaces divide the 
domain C K™ in subdomains Qj, with j = 1, • • • , k, of different region with specific internal parameter 
profiles. Due to the different physical structure of each of these regions, in general, quite different mathe- 
matical models might be most appropriate for describing them. Solutions of such model represents a free 
boundary problem, i.e., one in with interfaces are also unknown and must be determined in addition to the 
solution of the governing partial differential equation [9, 11, 8, 5, 10, 26, 41, 34], in general, by a set of data 
obtained by indirect measurements. Applications includes image segmentation problems [10, 26, 34, 41], 
optimal shape designer problems [32, 5], Stefan's type problems [5], inverse potential problems [15, 13, 14], 
inverse conductivity/resistivity problems [24, 17, 36, 11, 9, 5] among others [8, 11, 5, 32, 18]. 

There is often a large variety of priors information available for determining the unknown physical 
parameter, whose character depends on the given application. In this article, we are interesting in inverse 
problems that consists in the identification of an unknown quantity u € D(F) C X that represents all 
parameter profiles inside the individual subregions of £1, from data y £ Y, where X and Y are Banach 
spaces and D(F) will be adequately specified in Section 3. In this particular case, only the interfaces 
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between the different regions and, possibly, the unknown parameter values need to be reconstructed from 
the gathered data. 

This process can be formally described by the operator equation 



where F : D(F) C X — » Y is the forward operator. In practical applications, data are obtained only by 
indirect measurements of the parameter. Hence, in general, exact data y £Y are not known and we have 
only access to noise data y s 6 Y, whose level of noise 5 > are assumed be known a priori and satisfies 



In particular, problems as described above are ill-posed and the accuracy of an approximated solution 
call for a regularization method [19, 5]. In this article we propose a Tikhonov-type regularization method 
coupled with a level-set approach to obtain a stable approximation of the unknown level sets and values 
of the piecewise solution of (1). 

Many approaches, in particular level set type approaches, have been suggested for such problems 
previously. In [6, 7, 9, 22, 25, 30], level set approaches for identification of the unknown parameter u 
with distinct, but known, piecewise constant values was investigated. In [7, 10, 15], level set approaches 
was derived for solve inverse problems, assuming that u is defined by several distinct constant values. 
In both cases, one needs only to identify the level sets of u, i.e. the inverse problem reduces to a shape 
identification problem. On the other hand, when the level values of u are also unknown, the inverse 
problem becomes harder, since, we have to identify both, the level sets as well as the level values of the 
unknown parameter u. In this situation, the dimension of the parameter space increases by the number 
of unknown level values. Level set approaches to ill- posed problems with unknown constant level values 
appear before in [14, 13, 32, 33, 34]. For level set regularization properties of the approximated solution 
for inverse problems see [15, 13, 14, 22, 4]. However, regularization theory for inverse problems where 
the components of the parameter u are variable and has jumps of discontinuities was not well investigate. 
Some early results in the numerical implementation of such inverse problems was consider before in [11]. 

In this article, we propose a level set type regularization method to ill-posed problems whose solution 
is composed by piecewise components not necessarily constants. These type of problems, where the 
parameters are variable and discontinuous, arise in many problems as heat conduction between materials 
of different heat capacity and conductivity, interface diffusion processes and many other types of physical 
problems where modeling components are related with embedded boundaries. See, for example [8, 11, 18, 
3] and references. 

The advantage of our proposed level set method vis-a-vis the approaches before is that, we prove 
regularization properties for the approximated solution of the inverse problem model by (1) where the 
parameter is a non-constant piecewise solutions. Therefore such level set type method is different of those 
in the literature since we are able to apply this method to problems whose solutions depends of moving 
interfaces and non-constant parameters. In terms of applications, the presented regularization approach 
can handle with more kinds of flexibility and interesting physical problems, where de components of the 
interesting parameter is not necessarily homogeneous [18, 14, 13, 32, 33, 34, 3]. As a benchmark problem 
we analyze two inverse problem model by elliptic PDE's with discontinuous and variable interesting 
coefficient. 

The paper is outlined as follows: In Section 2, we formulate the Tikhonov functional based on the 
level-set framework. In Section 3, we present the general assumptions needed in this article and the 
definition of the set of admissible solutions. In Subsection 3.1, we prove relevant properties about the 
admissible set of solutions, in particular, convergence in suitable topologies. Relevant properties of the 
penalization functional is given in Subsection 3.2. In Section 4, we prove that the proposed method is a 
regularization method to inverse problem, i.e., we prove that the minimizers of the proposed Tikhonov 
functional are stable and convergent with respect to the noise level in the data. In Section 5, a smooth 
functional is proposed to approximate minimizers of the Tikhonov functional defined in the admissible 



F(u) = y, 



(1) 



y s - v\\y < s. 



(2) 
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set of solutions. We provide approximation properties and optimality condition for the minimizers of the 
smooth Tikhonov functional. In Section 6, we present an application of the proposed framework to solve 
some interesting inverse elliptic problem with variable coefficient. Conclusions and future directions are 
presented in Section 7. 

2 The Level-set Formulation 

Our starting point is the assumption that the parameter u in (1) assume two unknown functional values, 
i.e., u(x) G {ip 1 (x) , ip 2 (x)} a.e. in C R d , where fi is a bounded set. More specifically, we assume 
the existence of a mensurable set D CC with < \D\ < such that u(x) = ip l {x) if x G D and 
u(x) = ip 2 {x) if x G Vt/D. 

With this framework, the inverse problem that we are interesting in this article is the stable identifi- 
cation of both, the shape of D and the value function tjp(x) for x belongs to D and to fl/D, respectively, 
from observation of the data y s 6Y. 

We remark that, if ^(x) = c 1 and ip 2 (x) = c 2 with c 1 and c 2 unknown constants values, the problem 
of identify u was rigorously studies before in [14]. Moreover, many other approaches to this case appear 
in the literature. See [14, 7, 6, 5] and references. Recently, in [13], a L 2 level set approach to identify the 
level and constant contrast was investigated. 

Our approach differs from the level set methods proposed in [13, 14], by consider also the identification 
of variable unknown level of the parameter u. In this situation, many topological difficulties appears in 
order to have a tractable definition of an admissible set of parameters (see Definition 1 below) . General- 
ization to problems with more than tow levels following by the techniques derived in [15], applied to this 
approach. As observed before, the present level set approach is a rigorously derivation of a regularization 
strategies for identification of the shape and non-constant levels of discontinuous parameters. Therefore, 
it can be applied to physical problems modeled be embedded boundaries that whose components are not 
necessarily piecewise constant [22, 5, 15, 13, 14]. 

In many interesting applications, the equation (1) is ill-posed. Therefore, in order to obtain a stable 
approximate solution, some regularization method will be applied. We propose a regularization method 
by: First, introduce a parametrization on the parameter space, using a level set function <j> that belongs 
to H 1 ^). Note that, we can identify the distinct level sets of the function <p G H 1 ^) with the definition 
of the Heaviside projector 

H : if 1 ($7) — ► L°°(n) 

I other else . 

Now, from the framework introduced above, a solution of (1) can be represented as 

u{x) = ^ x {x)H{(j)) +V 2 (x)(l -H{<f>)) =: PO^^V 2 )^). (3) 

With this notation, we are able to determine the shapes of D as {x G 0; <f>(x) > 0} and Q/D as 
(f){x) < 0}. 

The functional level values ip 1 (x), ip 2 (x) are also assumed be unknown and they needed to be deter- 
mined as well. 

Assumption 1. We assume that ijj 1 ,^ 2 G IB := {/ : f is measurable andf(x) G [m, M] , a.e. inO}, for 
some constant values m,M. 

Remark 1. We remark that, f G IB implies that f G L°°(Q). Since O is bounded f G Moreover, 

J f(x)V ■ ip(x)dx < \M\ J \V-(<p)(x)\dx<\M\\\V-<p\\ Ll ( a yV<p€C%(fl,R n ). 
Hence f G BV(fi). 
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Note that, in the case that ip 1 and ip 2 assumes tow distinct constant values (as covered by the analysis 
done in [5, 14, 13] and references) the assumptions above are satisfied. Hence, the level set approach 
proposed here generalize the regularization theory developed in [14, 13]. 

From (3), the inverse problem in (1), with data given as in (2), can be abstractly written as the 
operator equation 

F(P(<^\Y> 2 )) = y S . (4) 

Since, an approximate solution (</>, ip , ip 2 ) of (4) is obtained, a corresponding solution of (1) can be 
computed using equation (3). 

Therefore, for obtain a regularized approximate solutions to (4), we shell consider the least square 
approach combined with a regularization term i.e., minimizing the Tikhonov functional 

G a {(p^\iP 2 ) := ||F(P(^,^\^ 2 ))-/||f, + a{/3 1 |^( ( /))|Bv + /32||^-</'o|lH 1 (n)+/ 5 3E I^'-^|bv}, (5) 

j=i 

where, cpo and ipQ represents some a priori information about the true solution u* of (1). The parameter 
a > play the rule of regularization parameter and the values of ,i = 1,2,3 act as scaling factors. In 
other words, ,i = 1,2,3 need to be chosen a priori, but independent of the noise level 5. In practical, 
fa , i = 1, 2, 3 can be chosen in order to represent a priori knowledge of features the of parameter solution 
u and/or to improve the numerical algorithm. A complete discussion about how to choose Pi ,i = 1,2,3 
are provided in [13, 14, 15]. 

The regularization strategy in this context is based on TV — H 1 — TV penalization. The term on H 1 - 
norm act simultaneous as a control on the size of the norm of the level set function and a regularization on 
the space H 1 . The term on BV is a variational measure of H((p). It is well known for penalizing the length 
of the Hausdorff measure of the boundary of the set {x : (p{x) > 0} (see [21, Chapter 5] for details). 
Finally, the last term on BV is a variational measure of ip 3 that acts as a regularization term on the set B. 
This Tikhonov functional extends the ones proposed in [14, 15, 6, 7, 32] (based on TV-H 1 penalization). 

Existence of minimizers for the functional (5), in the H 1 x B 2 topology do not follows by direct 
arguments, since, the operator P is not necessary continuous in this topology. Indeed, is tp 1 = tp 2 = tp a 
continuous function at the contact region, than P{cp 1 ,ip 2 ,tp) = ip is continuous and the standard Tikhonov 
regularization theory to inverse problem holds true [19]. On the other hand, in the interesting case where 
ip 1 and ip 2 represents the level of the jumps in the parameter u, the analysis became more complicated and 
we need a definition of generalized minimizers (see Definition 1) in order to handle with the difficulties. 

The advantage of our approach in relation to [5, 14, 22, 15] is that we are able to obtain regularize 
solutions to inverse problems with piecewise solutions that are more general that those covered by the 
regularization methods proposed before. In particular, we propose a level set regularization method to 
regularize piecewise non-constant solutions of ill-posed problems. Therefore, the approach studied here 
can be applied to a variate of physical problems of interesting as heat conduction between materials of 
different heat capacity and conductivity and interface diffusion processes, where the coefficient are variable 
and discontinuous [8, 11, 3]. 

We remark that, level set type methods was used for obtain solutions of elliptic problems with discon- 
tinuous and variable coefficients [11]. However, level set regularization theory for solve inverse problems 
[13, 15, 14] for recover the shape and the values of variable discontinuous coefficient are unknown be the 
author. 

The principal difference between our approach vis a vis the other ones before [14, 33, 32, 34, 30, 25, 
22, 9, 6, 7], is that we introduce a level set type regularization for recover the shape and the values of 
variable discontinuous coefficient. In this framework, the topologies needed to guaranties the existence 
of a minimizer (in a generalized sense) of the Tikhonov functional (5) is quite complicate and differ in 
some key points of [13, 14]. In these particular approach, the definition of generalized minimizer are quite 
different of others before [14, 22, 15]. As a consequence, the arguments used to prove the well-posed of the 
Tikhonov functional (5), the stability and convergence of the regularized solutions of the inverse problem 
(1) are quite complicated and needed significant improvements (see Definition 1). 
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3 Generalized Minimizers 

As already observed in [22], if D C ft with W 1 ' 1 ^) < oo where H n 1 (S) denotes the (n-l)-dimensional 
Hausdorff-measure of the set S, then the Heaviside operator H maps if 1 (ft) into the set 

V := { XD ; D C ft measurable, : U^idD) < oo} . 

Therefore, the operator P in (3) maps i7 1 (ft) x B 2 into the admissible parameter set 

D(F) := {u = q(y, i/j 1 , ip 2 ) ; v G V and ip 1 ,^ 2 

where 

g : V x B 2 9 (u,^ 1 ,^ 2 ) ^ ^^ + ^ 2 (1 - u) G W(ft) . 

Consider the model problem described as in the introduction. In this article, we assume that: 
(Al) ft C R n is bounded with piecewise C 1 boundary 5ft. 

(A2) The operator F : D(F) C L x (ft) — > Y is continuous on D{F) with respect to the L 1 (ft)- 
topology. 

(A3) e, a and ,j = 1,2, 3 denote positive parameters. 

(A4) Equation (1) has a solution, i.e. there exists G D(F) satisfying F(u*) = y and a function 
(/>* G -?/ 1 (ft) satisfying |V0*| 7^ 0, in a neighborhood of {0* = 0} such that H(<fi*) = z*, for some 
G V. Moreover, there exist functional values ipl,ip 2 G B such that q(z*, tpl , ip 2 ) = u*. 

For each e > 0, we define an smooth approximation to the operator P by 

P e (0,^,^) := ^ 1 F £ (0) + ^ 2 (1-F £ (0)) , (6) 
where -ff e is the smooth approximation to -ff describe by 



H e (t) :-- 



1 + t/e for t G [-e,0] 
£T(t) for t G E/[-e,0] ' 



Remark 2. is worth noting that, for any (pk G i/ 1 (ft), H e (<pk) belongs to L°°(ft) and satisfies < 
H £ (4>k) < 1 a.e. in ft, /or a// e > 0. Moreover, taking in account that ip 3 G B, follows that the operators 
q and P e , as above, are well defined. 

In order to guarantee existence of a minimizer of Q a defined in (5) in the space H l {pP) x B 2 , we need 
introduce a suitable topology such that the functional Q a has closed graphic. Therefore, the concept of 
generalized minimizers (compare with [15, 22]) in this paper is: 

Definition 1. Let the operators H , P, H £ and P £ be defined as above and positive parameters a,(3j and 
e given in the general assumption. 

A quadruple (z, (p, ip 1 , ip 2 ) G L°°(ft) x i^(ft) x BV(ft) 2 is called admissible when: 

a) There exists a sequence {4>k} of H 1 ^) -functions satisfying lim \\cpk — <^||,L 2 (r2) = 0- 

b) There exists a sequence {£k} G M + converging to zero such that lim \\H Ek ((f)j~) — -s||£i(n) = 0. 

c) There exist sequences {ip\}k&H an d {^i}keN belongs to BV n C°°(ft) such that 

I^Ibv^I^'Ibv, J = 1,2. 

d) A generalized minimizer ofQ a is considered to be any admissible quadruple (z, (p, ip\ ip 2 ) minimizing 

g a {z,^\^ 2 ) := \F{q{z,i>\^))-y s f Y + aR{zA^\^ 2 ) (7) 
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on the set of admissible quadruples. Here the functional R is defined by 



R(z,^ 1 ^ 2 ) = p(*,0)+ft£|^'-# 



2 

£ 

3=1 



and the functional p is defined as 



p(z, <f>) := inf { lim inf fa\H £k (4> k ) | BV + fa 
y. k— >oo 



OlBV 



9 I Iff l(O) 



(8) 



(9) 



The infimum in (9) is taken over all sequences {£k} and {4>k} characterizing (z, <f>, ij> , ip 2 ) as an admissible 
quadruple. 

The convergence |^|bv — ► IV^Ibv in Item c) in Definition 1 is in the sense of variation measure [21, 
Chapter 5]. The incorporation of item c) in the Definition 1 implies in the existence of the T-limit of 
sequences of admissible quadruples [22, 1]. This appears in the proof of Lemmas 4, 5 and 8, where we 
proof that the set of admissible quadruple are closed in the defined topology (see Lemmas 4 and 5) and 
the weak lower semi-continuity of the regularization functional R (see Lemma 8). The identification of 
non-constant level values ijj 3 implies in a different definition of admissible quadruple. As a consequence 
the arguments in the proof of regularization properties of the level set approach and the applications are 
the principal and novel difference between our definition of admissible quadruple in the present article 
and one that appears in [16, 13, 14, 22]. 



Remark 3. For j = 1, 2 let ^'eln C°°(f2), (f> E H 1 ^) be such that \V(f>\ ^ in a neighborhood of the 

t {(p(x) = 0} and H{<j>) = z ( 
°f { £ k}keN of positive numbers converging to zero, we have 



level set {4>(x) = 0} and H(<fi) = z S V. For each k € N set tpi = tpi and 4>k = 4>. Then, for all sequences 



\H, 



I i 1 (Q) 



\H, 



< 



e k •/(0)- 1 (r) 



\dr < meas{(</>)" 1 (r)} 



1 - — 



Idt 



dx 



Here, we use the fact that |V</>| ^ in a neighborhood of {<p = 0} implies that <j> is a local diffeomorphism 
together with a co-area formula [21, Chapter 4]- Moreover, {V'ilfceN inMnC°°(Q) satisfying Definition 1 
item c). 

Hence, (z, <f>, ip 1 , tp 2 ) is an admissible quadruple. In particular, we conclude from the general assump- 
tion above that the set of admissible quadruple satisfying F(u) = y is not empty. 



3.1 Relevant Properties of Admissible Quadruples 

Our first result will be to prove continuities properties of operators P e , H £ and q in suitable topologies. 
Such results will be necessary in the subsequent analysis. 

We start with a auxiliary lemma that is well known (see for example [16]). We present it here for sake 
of completeness. 

Lemma 2. Let Q be a measurable subset ofM n with finite measure. 

If (fk) SB is a convergent sequence in L P (Q) for some p, 1 < p < oo, then it is a convergent sequence 
in L p (f2) for all 1 < p < oo. 

In particular Lemma 2 holds for the sequence z^ := H e (cj)k). 

Proof. See [16, Lemma 2.1]. □ 

Next two lemmas are auxiliary results in order to understand the definition of the set of admissible 
quadruple. 
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Lemma 3. Let £1 as in the assumption (Al) and j = 1, 2. 

(i) Let {zk}keN be a sequence in L°°(Q,) with z k € [m,M] a.e. converging in the L l (Q)-norm to some 
element z and {V'fclfceN be a sequence in M converging in the BV-norm to some ip 3 £ B. Then 
q(zk,ipl,ip1) converges to q(z, ip 1 , ip 2 ) in L l (£l). 

(ii) Let (z,4>) £ L l (Q) x H l (£l), be such that H e (<j)) ->• z in L l (Vt) as e -> and let ip 1 ,ip 2 £ B. Then 
P^ip 1 ,^ 2 ) -> <?U V>\ V> 2 ) in L l (9.) as e ^ 0. 

(■m) Given e > 0, Ze£ {</>fc}fc<=N ^ e a sequence in H 1 ^) converging to £ H 1 ^) in the L?-norm. Then 
H £ (4>k) H e (<p) in L 1 (il), as A; — >• oo. Moreover, if {V'ilfceN are sequences in B, converging to 
some ip 3 in B, w'i/i respect to the L l (Q)-norm, then q(H £ ((pf c ), ipi, ipi) — > q(H e (cp), ip , ip 2 ) in L 1 (il), 
as k — ?> oo. 

Proof. Since is assumed to be bounded, we have L°°(r2) C L 1 (f2) and BV(£l) is continuous embedding 
in L 2 {Q) [21]. To prove (%), notice that 

- g(^,V' 1 ,V' 2 )llLi(n) = ||^fc«* + - z k ) - ip x z - ip 2 {l - z)\\ L1{n) 

- ll Z fcllL°°(n)ll^fc ~~ ^11 Li(H) + 11^ II L2(n)H Zfe ~ Z IIl 2 (C) 

.. ||2 211 II2II ^ 

+ W 1 - z k\\ L cv(n)\\iPk ~ "tp ILi(Q) + H"^ llL 2 (c)ll Zfc ~~ z \\l 2 (q) y °- 

Here we use Lemma 2 in order to guarantee the convergence of z k to z in L 2 (£2). 

Assertion (ii) follows with similar arguments and the fact that H e {<p) £ L°°{Q) for all e > 0. 

As \\H s ((pk) — H e ((p)\\ L i^ < meas(fi) \\cpk — 0||l 2 (Q) the first part of assertion (Hi) follows. The 
second part of the assertion ( Hi ) holds by a combination of the inequality above and steps in the prove of 
assertion (i). □ 

Lemma 4. Let {V'fclfceN be a sequence of functions satisfying the Definition 1 converging in L 1 (0) to 
some tp 3 , for j = 1,2. Then ip 3 also satisfies the Definition 1. 

Sketch of the proof. 

Let k £ N and j = 1,2. Since ip^, satisfies Definition 1, ipi £ BV. Prom [21, Theorem 2, pg 172] there 
exist sequences {ip J kl }i^ in BV x C°°(0) such that 

V'lZ-^'V'fc in L 1 ^) and |^|bv I^Ibv • 
In particular, for the subsequence {^ii^^keN follow that 

Vi/ (fe) ^</> j m L\n) and |^ (fc) | B v ^ IV^'Ibv • (10) 

Moreover, by assumption ip 3 £ L 1 (J7). From the lower semi-continuity of variational measure (see [21, 
Theorem 1 pg. 172]), equation (10) and definition of BV space follows that ip 3 £ BV. □ 
In the next lemma we prove that the set of admissible quadruple is closed with respect the L 1 (0) x 
L 2 (fl) x (L 1 ^)) 2 topology. 

Lemma 5. Let (z k , (pk^li^l) be a sequence of admissible quadruples converging in L X (J7) x L 2 (Q) x 
(L 1 (il)) 2 to some (z, <p, ip 1 , ip 2 ), with <p £ H 1 ^). Then, (z,<p,ip l ,ip 2 ) is also an admissible quadruple. 

Sketch of the proof. Let k £ N. Since {z\, (p\,ip\,ip\) is an admissible quadruple, it follows from Definition 1 
that there exist sequences {4>k,i}ieN, in H 1 ^), {ipl^i^, {ip 2 t }i^ in BV x C°°(f2) and a correspondent 
sequence {e^,}^ converging to zero such that 

<Pk,i l ^> 4>k in L 2 (Q), H £ i (cp k>l ) ^ z k in L l (£l) and |^|bv l= ^> IV^Ibv ,j = 1, 2 . 
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Define the monotone increasing function r : N — > N such that, for every k € N it holds 

e l {k) ^ l £ l~i 1) '. hk,T{k)-M\ L \a) ^ p ll^( fe )(^,-w)-^ILi(n) - p IV'v^Ibv — ► |V>fc|Bv,i = 1,2. 

(11) 

Hence, for each A; € N 

||0 - ^fe,r(fc)||i,2(f2) < 110 — 0fc||L2(n) + ||0fe, r (fc) _ 0fe||i2(Q) 
ll Z ~~ -^ e ^(fe)(0fc,r(fc))|| L l(Q) < \\Z - ZkWL 1 ^) + ||-^ e -r(fe)(0fe, r (fe)) — 2 fc|| L l(Q) • 



From (11) 



lim^ ||</> - 0ifc,r(fe)|| L 2(n) = , Hm o ||2-fl E r(fe)(0 fc)T ( fc ))||_ Ll( ^ =0. (12) 
Moreover, with the same arguments as Lemma 4, follows that 

l^fc,r(fe)l BV ~^ I^'IbV, J = 1,2, 

and S BV(J7). Therefore, it remains to prove that (z, 0, ip , tp 2 ) is an admissible quadruple. From 
Definition 1 and Lemma 4, it is enough to prove that z 6 L°°(f2). If this is not the case, there would exist 
a Q' C f£ with > and 7 > such that z(x) > 1 + 7 in (the other case: z(x) < —7 is analogous). 
Since (H T (k)(4>k T(k)))( x ) £ [0, 1] a.e. in for € N (see remark after Definition 1), we would have 

£ k 

\\z - H^(k)((f>k,T(k))\\Li(a) > lk--ff e r(fc)(0fc,r(fe))IUi(n') >7l fi 'l, fceN, 
contradicting the second limit in (12). □ 

3.2 Relevant Properties of the Penalization Functional 

In next lemmas, we verify properties of the functional R which are fundamental for the convergence 
analysis in Section 4. In particular, these properties implies that the level sets of Q a are compact in the 
set of admissible quadruple, i.e., Q a assume a minimizer on this set. First, we prove a lemma that simplify 
the functional R in (8). Here we present the sketch of the proof. For more details, see the arguments in 
[14, Lemma 3]. 

Lemma 6. Let (z, (f),^ 1 ,ip 2 ) be an admissible quadruple. Then, there exists sequences {£k}keN, {^fcjfceN 
and {ipi}keN as ^ n the Definition 1, such that 

RMrfrf) = hm I ^\H Ek (<j> k )\ BV + faWfa - M\m {Q ) + foJ2 W " ^oIbv J • (13) 

Sketch of the proof. For each / G N, the definition of R (see Definition 1) guaranties the existence of 
sequences e k , {4> J k z } £ H 1 ^) and {i(? k ( }el such that 



Now a similar extraction of subsequences as in Lemma 5 complete the proof. □ 
In the following, we prove tow lemmas that are essential to the proof of well posedness of the Tikhonov 
functional (5). 
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Lemma 7. The functional R in (8) is coercive on the set of admissible quadruples. In other words, given 
any admissible quadruple (z, <j), ip 1 , ip 2 ) we have 



Sketch of the proof. Let (z, cfrjip 1 ,tp 2 ) be an admissible quadruple. Prom [15, Lemma 4], it follows that 



p{z,<t>) > {Pi\z\w + foW - Mhpi)) ■ ( 14 ) 
Now, from (14) and the definition of R in (8), we have 

2 2 
3=1 3=1 



concluding the proof. 



□ 



Lemma 8. The functional R in (8) is weak lower semi- continuous on the set of admissible quadruples, 
i.e. given a sequence {(z^, ^|)} of admissible quadruples such that z k — > z in L 1 (il), — <f> in 
H l (Q), ip J k — > ip 3 in L 1 (fl), for some admissible quadruple (z, <p, ip 1 , ip 2 ), then 



R(z,(p,ip ,tp ) < liminf R(z k ,4>k,^ k ^k) ■ 

fcgN 

Proof. The functional p{z,<p) is weak lower semi-continuous cf. [15, Lemma 5]. As ip k £ BV follows from 



[21, Theorem 2 pg 172] that there exist sequences {ipi ; } G BV fl C°°(f2) such that ip kl — ip k 



From a diagonal argument, we can extract a subsequence {ip^n^} of {ip k i\ such that {ip^u^} — ► V'' 7 i n 
L 1 ^) as jfc -> oo. Let £ G C2(Q, M n ) , |f | < 1. Then, from [21, Theorem 1 pg 167], it follows that 

I' V • £dx = lim / ^ } V ■ £d* = lim / (fc) - ^) V ■ £dx + [ ^ V ■ £cfe 



< lim 



< liminf IV^Ibv 

fc— >oo 



Thus, form the definition of | • |bv (see [21]), we have 



W |bv = sup <^ / V J V • ; e G Q(ft, M n ) , |£| < l\ < liminf |^| BV . 

Now, the lemma follows from the fact that the functional R in (8) is a linear combination of lower 
semi-continuous functionals. □ 



4 Convergence Analysis 

In the following, we consider any positive parameter a, f3j ,j = 1, 2, 3 as in the general assumption to this 
article. First, we prove that the functional Q a in (7) is well posed. 

Theorem 9 (Well-Posedness). The functional Q a in (7) attains minimizers on the set of admissible 
quadruples. 
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Proof. Notice that, the set of admissible quadruples is not empty, since (0,0,0,0) is admissible. Let 
{(zk, 4>k, tpl, V'fc)} be a minimizing sequence for Q a , i.e. a sequence of admissible quadruples satisfying 
Q a {zk,4>k,ipk^k) ~~ ^ m f<?a — G Q (0, 0, 0, 0) < oo. Then, {Q a (zk, <Pk, Vi> ify} is a bounded sequence of 
real numbers. Therefore, {(z k , (f) k , ipi, V'fc)} is uniformly bounded in BV x x BV 2 . Thus, from the 

Sobolev Embedding Theorem [2, 21], we guarantee the existence of a subsequence (denoted again by 
{{z k ,<t> k ,4> k ,ip k )}) and the existence of (z, 4>, ip 1 , 4> 2 ) G L 1 (f2) x ff 1 (fi) x BV 2 such that fc ->• <^> in L 2 (f2), 
fc in ff 1 ^), z fc ^ z in L x (0) and ^ -> ^ J in L^fi). Moreover, z,^ 1 and V 2 € BV. See [21, 
Theorem 4, pp. 176]. 

From Lemma 5, we conclude that (z,<f>,ip ,ip 2 ) is an admissible quadruple. Moreover, from the weak 
lower semi-continuity of R (Lemma 8), together with the continuity of q (Lemma 3) and continuity of F 
(see the general assumption), we obtain 

inf g a = lim g a (z k ,4> k , ^1,^1) = lim {\\F(q(z k , - y S \\y + aR(z k , <j> k , tpl, tp k )} 

k— >oo k— >oo 

> WFiqiz^^-y'Wl + aRiz,^ 1 ^ 2 ) = Q a {z, <f>, ^, if) , (15) 

proving that (z, <j), ip , ip 2 ) minimizes Q a . □ 

In that follows, we shall denote a minimizer of Q a by (za^a^a-i^a)- In particular the functional Q a 
in (22) attain a generalized minimizer in the sense of Definition 1. In the next theorem, we summarize 
some convergence results for the regularized minimizers. These results are based on the existence of a 
generalized minimum norm solutions. 

Definition 2. An admissible quadruple (z' , <fv ,ip> 1 ^ ,ip 2,J( ) is called a R-minimizing solution if satisfies 
(i) F(q{z\^^)) = y, 

(ii) R{z\4>\^\i) 2 ^) = ms := inf {R(z,<f>, ip 1 , ip 2 ) ; (z,<j),ip l ,ip 2 ) is an 

admissible quadruple and F(q(z,ip 1 ,ip 2 )) = y} . 

Theorem 10 (R- minimizing solutions). Under the general assumptions of this paper there exists a R- 
minimizing solution. 

Proof. From the general assumption on this paper and Remark 3, we conclude that the set of admissible 
quadruple satisfying F(q(z, ip 1 ,ip 2 )) = y is not empty. Thus, ms in (ii) is finite and there exists a sequence 
{(z k , <pkt ^jfc> Y'/DlfeeN of admissible quadruple satisfying 

F {q( z k,i>l,ipl)) = y and R(z k ,4>k,^l,ipl) -> ms < oo . 

Now, form the definition of R, it follows that the sequences {</>fc}jfceNi {z k }keN and {V'fcligN 2 are uniformly 
bounded in H 1 (Q) and BV(f2), respectively. Then, from the Sobolev Compact Embedding Theorem [2, 21], 
we have (up to subsequences) that 

4>k^4> ] in L 2 (Q), z k ^z ] in L x (0) and ijP k ->• V i,f in L^O), j = l,2. 

Lemma 5 implies that (z\ <p\ , ip 2 ^) is an admissible quadruple. Since R is weakly lower semi- 
continuous (cf. Lemma 8), it follows 

ms = limmfR(zkAk^i^l) > R(z^ <f>\ ^ , ^ ) . 

k — ^oo 

Moreover, we conclude from Lemma 3 that 

q{z\^\^) = lim qMl^l) and F(q(z^,^\^)) = lim F(q(z k , ^)) = y . 

fc— S-oo k— >oo 

Thus, (z', <jy , ip 1 ^ i ip 2 '^) is a R- minimizing solution. □ 
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Using classical techniques from the analysis of Tikhonov regularization methods (see [20, 19]), we 
present below the main convergence and stability theorems of this paper. The arguments in the proof are 
somewhat different of that presented in [14, 13]. But, for sake of completeness, we present the proof. 

Theorem 11 (Convergence for exact data). Assume that we have exact data, i.e. y s = y. For 

every a > let (z a , cj) a ,ft,ft) denote a minimizer of Q a on the set of admissible quadruples. Then, for 
every sequence of positive numbers {a^keN converging to zero there exists a subsequence, denoted again 
by {afc}z e pj, such that (z ak ,<fi ak ,ft,ft) * s strongly convergent in L 1 (J7) x L 2 (S7) x (L 1 (Q)) 2 . Moreover, 
the limit is a solution of (1). 

Proof. Let (z*, 4> , ft', ft'*) be a i?-minimizing solution of (1) - its existence is guaranteed by Theorem 10. 
Let {afc}fcgN be a sequence of positive numbers converging to zero. For each k E N, denote (zk, 4>k, ft, ft) := 
{ z a k -, ( t , a k -,'^a k ->'^'a k ) be a minimizer of G ak . Then, for each k € N, we have 

G ak (z k , <j> k , ft, ft) < \\F(q(z\ft>\ ft'^)) - y\\ + a k R{z\ <j>\ft>\ ft A ) = a k R{z\ 4>\ ft*, ft*). (16) 

Since a k R(z k , , 4> k ,ft,ft) < G ak (z k ,(j) k ,ft,ft), it follows from (16) that 

R(z k ,<l> k ,ft,ft k ) < R(z\^,ft*,ft*) < oo. (17) 

Moreover, from the assumption on the sequence {a^}, it follows that 

lim a k R(z\ft,ft'\ft*) = 0. (18) 

k~ >OD 

From (17) and Lemma 7, we conclude that sequences {4>k}, { z k} and {ft k } are bounded in H 1 (Q)and 
BV, respectively, for j = 1,2. Using an argument of extraction of diagonal subsequences (see proof of 
Lemma 5), we can guarantee the existence of an admissible quadruple (z, <fr, ft, ft) such that 

(z k ,^ k ,ft,ft)^(z,4>,ft,ft) in L 1 (fl) x L 2 (S1) x (L 1 (il)) 2 . 

Now, from Lemma 3 (i), it follows that q(z, ft, ft) = lim q(z k ,ft,ft) in L 1 (0). Using the continuity of 

k— >oo 

the operator F together with (16) and (18), we conclude that 

y = lim F(q(z k ,ft,ft k )) = F(q(z,ft,ft)). 

k— >oo 

On the other hand, from the lower semi-continuity of R and (17) it follows that 

R(z,cf>,ft,ft) < ]immIR(z k ,(j) k ,ft,ft) < lim sup R(z k , <p k , ft, ft.)) < R(z\ <f>\ ft, ft) , 

k— ¥oo k-^oo 

concluding the proof. □ 

Theorem 12 (Stability). Let a = a(5) be a function satisfying lim ct(S) = and lim 5 2 a(5)~ 1 = 0. 

<5->o <5-»o 

Moreover, let {5fc}fceN be a sequence of positive numbers converging to zero and y Sk EY be corresponding 
noisy data satisfying (2). Then, there exist a subsequence, denoted again by {5 k }, and a sequence {a k := 
ct(S k )} such that {z ak ,4>a k ,ft k ,ft k ) converges in L 1 (il) x L 2 (Q.) x (L x (fi)) 2 to solution of (1). 

Proof. Let (z^, ft, ft*, ft*) beaR -minimizer solution of (1) (such existence is guaranteed by Theorem 10). 
For each k € N, let (z k ,4> k ,ft,ft) := (z a ( 5fc ), <t> a (s k ), ft a {s k y ^a{8 k )) be a mmimizer of G a(S k )- Then > for 
each k G N we have 

G ak (z k , <f> k , ft, ft) < \\F(q(zl ft*, ft*)) -y 5k \\ 2 Y + a(S k )R(z^,ft,ft*,ft*) 

< S 2 + a(5 k )R^,tf,ft*,ft*). (19) 
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From (19) and the definition of G ak , it follows that 

5 2 
a{Ok 



(20) 



Taking the limit as k — > oo in (20), it follows from theorem assumptions on a(5k), that 



lim 

k—>oo 



and 



F{q{z k ,^l))-y Sk < lim (5l + a(6 k )R(z\tf,^,^) 

k— >oo \ 



limsupi2(s fc ,0 fc ,Vjfc,^) < i?(^ t ,^ t ,^ 1 ' t ,V 2 ' t ) • 



0, 



(21) 



With the same arguments as in the proof of Theorem 11, we conclude that, at least a subsequence that 
we denote again by (zk,4>k,ipl,ip1), converge in L 1 (0) x L 2 (f2) x (L 1 (0)) 2 to some admissible quadruple 
(z, (ft, if} 1 , ip 2 ). Moreover, by taking the limit as k — > oo in (19), it follows from the assumption on F and 
Lemma 3 that 



F(q{z,4>,ip\ip 2 )) = lim F(q{z k ,ip\,ip 2 k )) = y. 



□ 



The functional Q a defined in (7) is not easy to handled numerically, i.e., we are not able to derive a 
suitable optimality condition to the minimizers of Q a . In the next section, we work in sight to surpass 
such difficulty. 



5 Numerical Solution 

In this section, we introduce a functional which can be handled numerically, and whose minimizers are 
'near' to the minimizers of Q a . Let G e ,a be the functional defined by 

QeM^rf) ■= \\F(P £ (<P,i}\if} 2 ))-y 5 \\ 2 Y + a{f3 1 \HM^ + f32n-Mm+^E I^'-^|bv), (22) 

3=1 

where P £ (<j>, ip ,ip 2 ) := q(H £ ((p),ip 1 , ip 2 ) is defined in (6). The functional Q E>a is well-posed as the following 
lemma shows: 

Lemma 13. Given positive constants a, e, /3j as in the general assumption of this article, <po G H (Q) 
and if>Q £ B, j = 1,2. Then, the functional Q e ^ a in (22) attains a minimizer on i? 1 (0) x (BV) 2 . 

Proof Since, inf{<? ei£t (0, ip 1 , if) 2 ) : (4>, ip 1 , ip 2 ) € H 1 ^) x (BV) 2 } < C? B , a (0,0,0) < oo, there exists a 
minimizing sequence {(4>k, ip k , ^1)} m H 1 ^) x B 2 satisfying 

lim QeAMWk) = hxi{g e , a {4>,iP\iP 2 ) : (W^^ff'^xB 2 }. 

k — ^oo 

Then, for fixed a > 0, the definition of Q £>a in (22) implies that the sequences {<pk} an d {^L} 3 2 
are bounded in i7 1 (0) and (BV) 2 , respectively. Therefore, from Banach-Alaoglu-Bourbaki Theorem [42] 
4>k — ^ 4> m H 1 ^) and from [21, Theorem 4 pg. 176], ipi — > ip 3 ' in L 1 (J7), j = 1, 2. Now, a similar argument 
as in Lemma 4 implies that ip 3 G B, for j = 1,2. Moreover, by the weak lower semi-continuity of the 
i^-norm [42] and | ■ | B v measure (see [21, Theorem 1 pg. 172]), it follows that 

110 - 00 1 1 h i < hminf \\(f> k - 0o||#i and \ip j - V>q|bv < liminf \ipi - V'qIbv ■ 

k—^oo k— >oo 

The compact embedding of H 1 ^) into L 2 (Q) [2] implies in the existence of a subsequence of {4> k }, 
(that we denote with the same index) such that <p k — > <p in L 2 (f2). Follows from Lemma 3 and [21, 
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Theorem 1, pg 172] that \H £ ((/))\ B y < liminf \H £ ((/) k )\ B y. Hence, from continuity of F in L , continuity of 

fc— ¥00 

q (see Lemma 3), together with the estimates above, we conclude that 

GeAM 1 ^) < , lim \\F(Pe(Ml,i>l))-y s \\Y 



k— ¥00 



+ a (fii liminf \H £ {4> k )\ m + j3 2 liminf \\4> k - (po\\ 2 H i(n) + £ 3 liminf Y, \^k ~ V>oIbv) 

V k— >oo k— too v ' k~ too j=\ ' 

< liminf Ge^k^l^D = ini G £ , a , 

k— too 

Therefore, ((p,ip ,tp 2 ) is a minimizer of G £ , a - d 

In the sequel, we prove that, when e — >• 0, the minimizers of Ge,a approximate a minimizer of the 
functional Go,- Hence, numerically, the minimizer of Gea can be used as a suitable approximation for the 
minimizers of Go- 

Theorem 14. Let a and /L be given as in the general assumption of this article. For each e > 0, 
denote by (0 £)Q; , a , ip 2 a ) a minimizer of G e ,a (that there exist form Lemma 13). Then, there exists 
a sequence of positive numbers e k — > such that (H £k (4> £kt a), 4> £k ,a-, i } \ k , a -> ^e fc ,a) converges strongly in 
L 1 (Q) x L 2 (f2) x (L 1 (17)) 2 and the limit minimizes G a on the set of admissible quadruples. 

Proof. Let (z a , cf) a , ip^, tp^) be a minimizer of the functional G a on the set of admissible quadruples (cf. 
Theorem 9). From Definition 1, there exists a sequence {e k } of positive numbers converging to zero 
and corresponding sequences {4> k } m H 1 ^) satisfying <p k — > 4>a in L 2 (ft), H £k (cp k ) — > z a in and, 
finally, sequences {ip k } in BV x C^°($l) such that |Y>{|bv — * IV^Ibv • Moreover, we can further assume (see 
Lemma 6) that 



Riza^^ip^ipl) = lim (Pi\H ek ((p k )\ BV + P2\\4>k ~ (poWmtn) + /%E Wk ~ V'oIbv) • 

k— too y ' 7=1 



Let ((p £k ,ip\ k ,ip 2 k ) be a minimizer of G £k , a - Hence, (4> £k ,ip\ k ,ip 2 k ) belongs to i^ 1 (0) xl 2 (see Lemma 13). 
The sequences {H £k (cp £k )},{(j) £k } and {ipi k } are uniformly bounded in BV(fi), H 1 ^) andBV(fi), for j = 1,2, 
respectively. Form compact embedding (see Theorems [2] and [21, Theorem 4 pg. 176]), there exist 
convergent subsequences whose limits are denoted by z, <p and ipfl belong to BV(il),// 1 (J7) and BV(f2), for 
j = 1,2, respectively. 

Summarizing, we have (p £k — > <p m -^ 2 (^)> H £k ((p £k ) — > z in L 1 (f2), and ipi k — >■ ip J in L 1 (r2), j = 1,2. 
Thus, € L\n) x H^Q) x L 1 ^) is an admissible quadruple (cf. Lemma 5). 

From the definition of R, Lemma 3 and the continuity of F, it follows that 

\\F{q{z^\^))-y s \\l = \^ k ^\\F(P £k ^ £k ^l k ^ 2 £k ))-y 5 \\ 2 Y , 

R(zJ,^\^ 2 ) < liminf fe ^ oo (/3i|iI efe (0 £ J|Bv + /32||</> £fc -0o|||i (Q) +/33E|=il^ fe -'^|Bv) • 

Therefore, 

Ga(zA,^,i> 2 ) = \\F(q(z,^\^ 2 ))-y S \\^ + aR(z,^ 1 ,i> 2 ) 

< liminf G £kta (<j) £k ,ip 1 ip 2 ) < liminfg £ ^(fa^l^l) 

k-too K K k— too 

< hmsup||F(P £fc (^,^,^))-/l|y 

k— too 

+ alimsup (j3i\H ek ((p k )\ BV + fcWfa ~ ^o\\h^{U) + &£f=i Wh ~ V'oIbv) 

k— too 

= \\F{q{z a ,^i^l))-y S \\ 2 , + aR{z a ,(p OL ^l L ,'ip 2 a ) = G a (z a , <t>l, tpl, tpl) , 
characterizing (z, (p,ipa,ipa) as a minimizer of Go- D 
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5.1 Optimality Conditions for the Stabilized Functional 

For numerical purposes it is convenient to derive first order optimality conditions for minimizers of the 
functional Q a . Since P is a discontinuous operator, it is not possible. However, thanks to the Theorem 12, 
the minimizers of the stabilized functionals Q e ^ a can be used for approximate minimizers of the functional 
Q a . Therefore, we consider Q £a in (22) ,with Y a Hilbert space, and we look for the Gateaux directional 
derivatives with respect to eft and the unknown ip 3 for j = 1,2. 

Since H' e ((j)) is self-adjoint 1 , we can write the optimality conditions for the functional Q £ a in the form 
of the system 



a(A-I)(<f>-<f>o) = L £>a ^((f), ip 1 ,^ 2 ) , in ft (23a) 

(0-0o) -V = 0, at dn (23b) 

a y [V(^-^o)/|V(^'-^)|] = 4^,^,^), i = 1 > 2 - ( 23c ) 
Here v{x) represents the external unit normal quadruple at x G d£l, and 

-/3 1 (2/3 2 )- 1 ^(^) V • [Vtf e (<«/|Vtf e (0)|] , (24a) 

Ll^i^ 1 ^ 2 ) = (2/3 3 )" 1 (F / (P e (^,^ 1 ,V 2 ))^(^)*(i ? (P(^V' 1 ^ 2 ))-/) (24b) 
Ll a3 {^\^) = (2/33)- 1 (F / (P e (^^ 1 ,V 2 ))(l-^(^)))*(P(P(^V' 1 ,V' 2 ))-/). (24c) 



It is worth noticing that the derivation of (23) is purely formal, since the BV seminoma is not differen- 
tiable. Moreover the terms \VH e (^>)\ and \S/{ip 3 — iPq)\ appearing in the denominators of (23) and (24), 
respectively. 

In Section 6, system (23) and (24) is used as starting point for the derivation of a level set type 
method. 

6 Inverse Elliptic Problems 

In this section, we discuss the proposed level set approach and their application in some physical problems 
model by elliptic PDE's. We also discuss briefly the numerical implementations of the iterative method 
based on the level set approach. We remark that the iterative algorithm derived by the level set approach 
need a early stoop criteria [19], in the case of noise data. 

6.1 The Inverse Potential Problem 

In this subsection, we apply the level set regularization framework developed before in this article in an 
inverse potential problem, similar to ones considered in [38, 39, 13, 14, 22, 36, 23, 37]. 

The forward problem consists of solving on a given Lipschitz domain Q C M. n , for a given source 
function u G L 2 (Q) and a boundary function g G L 2 (dQ), the Poisson boundary value problem 

- V • (ffVw) = u, in ft , jiw + 72ii>i/ = g on ffi. (25) 

In (25), v represent the outer normal vector to d£l, a is a known sufficient smooth function. Note that, 
depending of 71, 72 G {0, 1}, we have Dirichlet, Neumann or Robin boundary condition. In this paper, we 
only consider the case of Dirichlet boundary condition, that corresponds to 71 = 1 and 72 = in (25). It 
is well known that there exists a unique solution w G i? 1 (0) of (25) with w — g G Hq(Q), [12]. 

These simple version of the elliptic problem (25) appears in many relevant applications including: 
inverse potential problems [23], inverse gravimetry [24] and EMG [39, 38]. 



1 Note that H' e (t) = K ' ' 

other else . 
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Assuming homogeneous Dirichlet boundary condition in (25), the problem can be modeled by the 
operator equation 

F-y : L 2 {n) -4 L 2 (dn) 

u\ — > F\ (u) := w„\on . (26) 

The corresponding inverse problem is, in general, called inverse potential problem, which consists of 
recovering an L 2 -function u, from measurements of the Cauchy data of its corresponding potential on the 
boundary of Q. 

Using this notation, the inverse potential problem can be written in the abbreviated form F\(u) = y s , 
where the available noisy data y s £ L 2 (dfl) have the same meaning as in (2). It is worth noticing that 
this inverse problem has, in general, non unique solution [23]. Sufficient conditions for identifiability are 
given in [24]. Moreover, in this case, the operator F\ is linear. However, the inverse potential problem 
is well known to be exponentially ill-posed [24]. Therefore, the solution call for a regularization strategy 
[19, 23, 24]. 

Given the non-uniqueness of the inverse potential problem and the limitation of our approach, we 
restrict our attention to solve the inverse problem (26) in D(F), i.e., we assume that the unknown 
parameter u £ D(F), as defined in Section 3. For this class of parameters no unique identifiability result 
is known and we restrict our attention to minimum-norm solutions [19]. The following lemma implies 
that the operator Fi satisfies the Assumption (A2). 

Lemma 15. The operator F\ : D(F) C L 1 (f2) — > L 2 (dQ) is continuous with the respect to the L 1 (J7) 
topology. 

Proof. It is well known form the elliptic regularity theory [12] that IMI.f/i(Q) ^ c ilMlL 2 (n)- Let u m u o £ 
D(F) and w n ,wo the respective solution of (25). Then, the linearity and continuity of the trace operator 
from H 1 ^) to L 2 (dVt) [12], we have that 

\\Fi(u n ) - Fi(n )|| ^(an) ^ C\\w n - w \\ H i^ < C\\u n - u \\ L2 ^ . 

From Lemma 2, we conclude that F\ is sequentially continuous on the L l (Q) topology. Since L l (H) 
is a metrizable spaces [42], the proof is complete. □ 

6.1.1 A level set algorithm for the inverse potential problem 

The iterative algorithm based on the level set approach proposed in this article is an explicit iterative 
method derived from the optimality conditions for the Tikhonov functional Q Ei a m (23) and (24). To the 
inverse potential problem for Dirichlet boundary condition (71 = 1 and 72 = 0) the algorithm reads as: 

Each step of this iterative method consists of three parts (see Table 1): 1) The residual € L 2 (d£l) 
of the iterate V^) is evaluated (this requires solving one elliptic BVP of Dirichlet type); 2) The L 2 - 
solution hk of the adjoint problem for the residual is evaluated (this corresponds to solving one elliptic 
BVP of Dirichlet type); 3) The update 54>k for the level-set function and the updates 5ip 3 k for the level 
values are evaluated (this corresponds to multiplying two functions). 

In [37], a level set method was proposed, where the iteration is based on an inexact Newton type 
method with the inner iteration is implemented using the conjugate gradient method. Moreover, the 
regularization parameter a > is kept fixed. In the above algorithm, we define 5t = 1/a as time 
increment in order to derive an evolution equation for the levelset function. Therefore, we are looking for 
a fixed point equation related to the system of optimality conditions for the Tikhonov functional. Here 
the iteration is based on a gradient type method as in [13]. 
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Given a and g; 

1. Evaluate the residual r k := Fi(P £ ((p k , - y 5 = {w k ) v \dn ~ y 5 > where w k solves 

-V • (crVw fc ) = P E {(j) k ,i)l,il>l) , in Q; w k = g , at dU . 

2. Evaluate h k := F{{P e (<j> k , ^, e . solving 

A/ifc = , in ; h k = r k , at d£l . 

3. Calculate b$ k \= L e ^p{§ k ^\^ 2 k ) and <^ := L 3 ^ aj3 (<j) k , ^) , as in (24). 

4. Update the level set function 0/% and the level values ^i, j = 1,2: 

<t>k+\ = «fo + £ , tf k+1 =ji + hML 

Table 1: Iterative algorithm based on the level set approach for the inverse potential problem. 
6.2 The Inverse Problem in Nonlinear Electromagnetism 

Many interesting physical problems are model by quasi-linear elliptic equations. One such problem that 
is interesting in applications is the identification of inhomogeneity inside nonlinear magnetic materials 
form indirect or local measurements. Electromagnetic non-destructive tests aim to localize cracks or 
inhomogeneities in the steel production process where the impurities can be described by a piecewise 
smooth function, [8, 9, 5, 11]. 

In this section, we assume that D CC 0, is measurable and 

)ipi,x € D , 
u = < 

[^2,x Gn D , 

with ipi,ifj2 S IB with m > 0. 

The forward problem consists of solving on a given Lipschitz domain Q C M n , for a given source 
function / € a boundary function g £ H l l 2 (dQ), the Poisson boundary value problem 

- V • [uVw) = f , in n , w = g on dfl. (28) 

It is well known that there exists a unique solution w € H l (£l) such that w — g G Hq(£2) for the problem 
(28), [12]. 

Assuming that, during the production process, the workpiece is contaminated by impurities and that 
such impurities are described by piecewise smooth function, the inverse electromagnetic problem consist 
in the identification and the localization of the inhomogeneities as well as the function values of the 
impurities. The localization of the support of the inhomogeneities and the tabulation of his values can 
indicate possible sources of contamination in the magnetic material. Note that, if ipi and ip2 are constant, 
then (28) describe the EIT problem [36, 24, 6, 7]. 

The inverse problem that we are interesting here consists in the identification of the support (shape) 
and the function values of tp 1 , i/j 2 of the coefficient function u(x) defined in (27), where the voltage potential 
g is chosen such that its corresponding the current measurement h := (w) u \gn are available as a set of 
continuous measurement in d£l. This problem is known in the literature as the inverse problem for the 
Dirichlet-to-Neumann operator [24]. 

With this framework, the problem can be model by the operator equation 

F 2 : D(F) C L 1 ^) -> H 1/2 (dn) 

u .— > F 2 (u) := w\ dn , (29) 

where the potential profile g = w\gQ € H l l 2 {Q) is given. 
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In [11], a level set approach for solve an inverse problems of identification of inhomogeneities inside a 
nonlinear material form local measurements of the magnetic induction was investigate. The assumption in 
[11] is that part of the inhomogeneities are given by a crack localized inside the workpiece and that outside 
the crack region, magnetic conductivities are nonlinear and depends on the magnetic induction. In other 
words, that ip\ = \i\ and ip 2 = A*2(|Vu;| 2 ), where /ii is the (constant) air conductivity and \i 2 = ^2(|Vu>| 2 ) 
is a nonlinear conductivity of the workpiece material, whose values are assumed be known. A successful 
iterative algorithm and numerical experiment was presented. However, we remark that, in [11], the 
measurements and therefore the data are given in the whole 0. Such measurements are not reasonable 
in applications. Moreover, the proposed level set algorithm is based on an optimality condition of a least 
square functional with H 1 (ri)-semi-norm regularization. However, there is no guarantee of existence of 
minimum for the proposed functional. 

In the following lemma, we prove that the operator F 2 satisfies the Assumption (A2). 

Remark 4. Note that the operator F 2 (u) = Td(w), where Tq is the Dirichlet trace operator. Moreover, 
Td : H 1 ^) — > iT 1 / 2 ^^) is linear and continuous [12], i.e., \\Td (w) Wh 1 / 2 ^) — c IIHI.H"i(n)- 

Lemma 16. Let the operator F 2 : D(F) C LP(9,) — > H 1 / 2 ^) as defined in (29). Then, F 2 is 
continuous with the respect to the L 1 (f2) topology. 

Proof. Let u n ,uo € D(F) and w n ,wo denoting the respective solution of (25). The linearity of equation 
(28) implies that w n — wq € Hq(Q) satisfies 

V • (u n Vw n ) - V • (u Vw ) = , (30) 

with homogeneous boundary condition. Therefore, using the weak formulation for (30) we have 

/ (V • (u n Vw n ) - V • (u Vw )) ipdx = , V</? € Hq(TI) . 
Jn 

In particular, the weak formulation holds true for ip = w n — wq. From the Green formula [12] and the 
assumption that m > (that guarantee elipticity of (28)), follows that 

m\\Vw n -\7w \\ 2 L 2 (n) < / u n \Vw n - Vw \ 2 dx < / \(u n - u )\\Vw \\(Vw n - Vw )\dx . (31) 
y ' Jn Jn 

From [28, Theorem 1], there exist e > (small enough) such that wo G W 1,P (Q) for p = 2 + e. Using the 
Holder inequality [12] with l/p+ 1/q = 1/2 (note that q > 2 in the equation (31), follows that 

m\\Vw n -Vw \\ 2 L 2 iu) < \\un-u \\ Lq{n) \\Vw \\ LP{n) \\Vw n -Vw \\ LHn) . (32) 

Therefore, using the Poincare inequality [12] and equation (32), we have 

\\ w n - w o\\ H i(n) - C\\ u n - ""oil L 9(fi) , 

where the constant C depends only of m, f2, ||ViOo[| an d the Poincare constant. Now, the assertion follows 
from Lemma 2 and Remark 4. □ 

6.2.1 A level set algorithm for inverse problem in nonlinear electromagnetism 

The explicit iterative algorithm based on the level set approach proposed in this article is based on the 
optimality conditions for the Tikhonov functional Q ea in (23) and (24). Each iteration of this algorithm 
consists in the following steps: In the first step the residual vector r £ L 2 (dCl) corresponding to the iterate 
(4> n , 11)^, ip 2 ) is evaluated. This requires the solution of one elliptic BVP's of Dirichlet type. In the second 
step the solutions v € H 1 ^) of the adjoint problems for the residual components r are evaluated. This 
corresponds to solving one elliptic BVP of Neumann type and to computing the inner-product Vw ■ Vi> in 
L 2 (f2). Next, the computation of L £t a t p((j) n , ip^, ip^) an d L 3 e a n[<p n ,ip\, ip^) as in (24). The four step is the 

updates of the level-set function 5cp n € H 1 ^) and the level function values Sipi, € BV(f2) by solve (23). 
The algorithm is summarized in Table 2. 
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1. 


Evaluate the residual r := F2{P E (4>m V^j V^)) ~ V & = w \dn ~ 5 ,<5 > where id 
solves 

V-(P £ (<^,^,^) V w) =/, infi; = at 90. 


£ tf 1 ^) 


2. 


Evaluate F,£(P e (</> n , ^))* r := Vw ■ Vv £ L 2 (0), where w is the 
computed in Step 1. and v £ i? 1 (0) solves 

V-(P e (0n,^,^)Vu) =0 in 0; ^ = r , at dtl . 


function 


3. 


Calculate L e ^ fi {^> n ,^\,^ 1 n ) and L J £ ^ (<p n , ij>*, ip%) as in (24). 




4. 


Evaluate the updates o<p £ il (S2) , of/v £ BV(O) by solving (23) 




5. 


Update the level set functions </> n +i = n + — 8<f), and the level 
values ^n+i =4> J n+ ■ 


function 



Table 2: An explicit algorithm based on the proposed level set iterative regularization method. 



7 Conclusions and Future Directions 

In this article, we generalize the results of convergence and stability of the level set regularization approach 
proposed in [14, 13] to the case where the level values and the level set are non-constant inside of each 
region. We analyze the particular case, where the set is divide in two regions. It is easy to extend the 
analysis to the case of multiple regions adapting the multiple level set approach in [15, 14]. 

We shows the applicability of the level set framework introduced in this paper for an inverse potential 
problem with piecewise smooth components. The application is a natural generalization of the inverse 
potential problem computed in [15, 13, 14]. Applications of inverse potential problems in EMG appears 
in [39, 38]. We also investigate the applicability in the identification of inhomogeneities inside a nonlinear 
magnetic workpiece. In both application, an iterative algorithm based on the optimality condition of the 
smooth Tikhonov functional Q EtCC was proposed. 

A natural continuation of this paper is the numerical implementation. Level set numerical imple- 
mentations for the inverse potential problem was done before in [15, 14, 13], where the level values are 
assumed to be constant. Implementations of level set methods for resistivity/conductivity problem in 
elliptic equation have been intensively implemented recently. See, for example, [9, 18, 31, 40, 36, 11, 5] 
and references. 
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